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Abstract 

We prove known identities for the Khinchin constant and develop new 
identities for the more general Holder mean limits of continued fractions. 

Any of these constants can be developed as a rapidly converging series 
involving values of the Riemann zeta function and rational coefficients. 

Such identities allow for efficient numerical evaluation of the relevant con- 
stants. We present free-parameter, optimizable versions of the identities, 
and report numerical results. 
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1 Introduction 


The Khinchin constant arises in the measure theory of continued fractions. Ev- 
ery positive irrational number can be written uniquely as a simple continued 
fraction [a 0 , aj, 02 , • • • , a„, • • •]; i.e., with a 0 a non-negative, and all other a, pos- 
itive integers. The Gauss-Kuz’min distribution ([11]) predicts that the density 
of occurrence of some chosen positive integer k in a random such fraction is 
given by 

Prob(a n = k) = - log 2 Jl - — ^ 2 . 

In his celebrated text, Khinchin ([11]) uses the Gauss-Kuz’min distribution to 
show that for almost all positive irrationals the limiting geometric mean of the 
positive elements a,* of the relevant continued fraction exists and equals 


A'. := n 


k — l 


1 + 


1 

k(k + 2). 


log 2 k 


JJ jfc lo e4 1+ TnW] . 

k = 1 


The fundamental constant Ao is the Khinchin constant . It is known that this 
constant can be cast in terms of various converging series, the following exam- 
ple of which having been used decades ago to provide the first high-precision 
numerical values for K 0 [16, 21, 22] 


log(A'o)log(2)^f; C(25) 1 

,=i s 



2s - 1 


( 1 ) 


This series can be rendered even more computationally efficient via the introduc- 
tion of a free integer parameter. We used a carefully optimized free- parameter 
series to resolve K 0 to over 7000 decimal places (A' 0 = 2.68545200106 •• • see 
Section 5). 


The Khinchin constant can be thought of as a special case of a more general 
statistical mean. For any real number p < 1, the Holder mean of order p of the 
continued fraction elements, namely lim t [(a^ + a p 2 + ■ ■ ■ + a*)/fc] 1/p , also exists 
with probability one and equals 


h p 


OO 

Y! ~ hP lo S2 


1 


1 

(Jfe + 1) 2 . 


1/p 


(See the final section of Khinchin *s book ([11]) for a proof for p < or more 
modern references on ergodic theory for a proof for p < 1 [15].) We may interpret 
I<o as the limiting instance of the K p definition as p — ► 0. We shall show that 
for any negative integer p the Khinchin mean of order p satisfies an identity 


(R p )Mog(2) = £(C(s-p)-l)Q Jp 

5 = 2 


( 2 ) 
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where each coefficient Q sp is rational. Again there is a free-parameter general- 
ization, which we ein ployed to resolve the harmonic mean A'_i also to over 7000 
decimal places {K\ — 1.74540566240- •• see Section 5). It is of interest that, 
evidently, only A'o can be written as a series involving exclusively even zeta 
arguments. The computational implications of this unique property of A'o are 
discussed in Section 5. We should mention that aside from the Shanks- Wrench 
series for A'o there are other previously known formulae for Khinchin means, 
some of which formulae involving derivatives of the zeta function [19]. 


2 Fundamental Identities 


This section is devoted to presenting the basic identities. We begin with a list 
of preliminary, largely elementary, results needed in the paper. All of these 
rearrangements may be easily justified. 

Lemma 1. (a) 

o° - 

- log( 1 - X) log(l + *) = ^2 -J-X 2k 

k = 1 

where A s := 

(b) 


Elogd - ^)log(l + 1) Tlog(* - l)log(l - ~ \ = -log(A)log(l + ~). 

k = 2 k—2 

Thus, (c) 

°° 1 1 

£log(l - j;) lo g(l + = -log(A'o)log(2). 

k = 2 

Proof. Part (a) is most easily seen by differentiating both sides. The left-hand 
side becomes f(x) — /(— x) where /(x) := log(l +x)/(l — x). Using the standard 
relationship 


oo 



k = 1 


oo k 




x 


k 


produces (a). 

Part(b) is easily established inductively after expanding the left-hand side. 
Part (c) follows on taking limits and noting that 


N 

~E l ° g(* _ FMi 

k=:> 


£ 2 ) = log(A' 0 )log(2), 
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as follows from the definition of A' 0 . 


0 


We will find it convenient to use the Hurwitz zeta function written 


C(*,N) 


= E 


1 

(n + N) s 


so that £(s) — £(s, 0) and so that for N a non-negative integer 


«*, ")=«•) -it ^ 

n = 1 

With this notation we have: 

Lemma 2. (a) For N a positive integer 


n= 1 


1 

~N' 


(b) For N a positive integer 


E 


<(2 n, AQ 
n 


= log(N + 1). 


(c) 



log(l - t 2 ) 
t( 1 -b 2) 


dt = - log 2 (2). 


Proof. The proofs of the first two identities are similar and rely on expanding 
the zeta terms, rearranging the order of summation and re-evaluating. In both 
cases, the result telescopes to the desired conclusion. 

Part (c) is less immediate. Actually, the indefinite integral is evaluable with the 
aid of the dilogarithm ([12]). Alternatively, Maple yields the result that 


equals the log terms 


Jo *(!+*) 


log 2 (l + „ 

2 log ( 2 ) + lo S( 2 ) lo S( l ~ t )~ i°g( 1 + 0 log(l - 1) + log(i) log( 1 - 1) 

plus the dilog terms 


dilog(0 - dilog(l + t) - dilog(iy^), 

and Maple also happily performs the final evaluation (the limit at 1). (L) 
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We are now in a position to establish the general Shanks- Wrench identity [16] 
for Ko- 

Theorem 3. For any positive integer N , 


oo . N . 1 

log(A'o)log(2) = - J>g(l - 7)log(l + T). (3) 

s = 1 S k-2 

where A t := ZVj'i -l)^ 1 /*• 

Remark. jV is a free parameter that can be optimized in actual computations 
to significantly reduce the number of zeta evaluations required. Variation of 
this parameter also provides a kind of error check, for whatever the choice of 
positive integer N , one expects an invariant result for the left-hand side. Note 
that in the case A r = 1 the second summation is empty, and we recover precisely 
the A'o identity (1) of Section 1. 

Proof. Let f(N) denote the right-hand side of (1). Then 

o° A 1 , 

f(N - 1) - /(AO = £ + log(l - -) log(l + -) 

5 = 1 

which equals zero by Lemma 1(a). Thus, since £(2s, N) — ► 0, sufficiently rapidly 


/(l)=/(yv)= lim f(N) 

N-^oo 

00 i | 

= ~IZ lo g(i - r)iog(i + t)- 

k = 2 

By Lemma 1(c), this sum agrees with log(A'o) log(2). (C) 

As a companion relation to the identity of Theorem 3, we can establish an 
elegant integral representation for the left-hand side. There is a powerful gen- 
eralization of Lemma 2(b) in the form of a generating function based on Euler’s 
product for s\n(7rt)/(irt) (see [18], page 249). For real t in [0,1) define g(t) by 


»(<) : = E 


C(2s) — 1 


t 2s = - log 


“"<">) -HogO-,’). 


7 rt 


(4) 


and define also the limiting case g{\) ;= log 2. We only need observe now that, 
on the basis of Theorem 3, with parameter N = 1, 


>g(/to)log(2) = [ 
Jo 


l°g(2) + g(t)/t 

1 + t 


dt. 


and with the help of the previous dilogarithm integral evaluation we thus arrive 
at an integral representation. ([16] contains an equivalent, though less stream- 
lined, integral identity.) 
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Corollary 4. The following integral representation holds for K 0 : 


l 


1 log[sin(7r<)/ (rtj\ 

^(1 + t ) 


dt = - log(A'o) log(2). 


The integral is sufficiently easy to handle that Greg Fee has computed 500 digits 
of A' 0 in Maple this way in about 10 minutes. It is amusing to observe that 
Lemma 1(c) may also be turned into an an analogous integral form: 


log(A' 0 )log(2)= 1“ = f 

Ji n 1 + 1) j 0 i + 1 

This was observed from a very different starting point by Robert Corless [9] but 
follows immediately on breaking the first integral up at integer points, 

We now derive new, corresponding identities for the higher-order Khinchin 
means. They are in some sense simpler, since one logarithmic term is replaced 
by a negative integral power. There is an observation that leads directly to a 
zeta function expansion for these general Khinchin means. Note that a sum of 
terms k p log(l — (k + l)” 2 ) can be expressed, via expansion of the logarithm, in 
terms of sums of the form (note p is assumed to be a negative integer): 


E 


1 

n 2 *-P(l - 1 /n)~P ‘ 


Upon expansion of the term 


1/(1 -1/n)-' 


in powers of 1/n, we obtain an identity for the p - th power of K p as a series 
of zeta functions. The result, after the same free parameter manipulations we 
used for K 0y reads: 

Theorem 6. For negative integer p and positive integer N we have 


Kp Iog(2) = £ 

n = l 


Er=o + 


N 


E 1 ^ 1 "pH* -0 P 


k — 2 


Remark. Note that for N — 1 the final sum is empty, the coefficient of any 
given C(s) is an easily computed rational, and we immediately establish a general 
series with rational coefficients, (2) of Section 1. 

Corollary 7. The harmonic Khinchin constant satisfies for integer N > 0 : 

l°g(2) _ N ~ N) _ log( 1 - &~ 2 ) 

K ~ i At n h k ~ i ' 
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Proof. It suffices to show that 

o 2 n . 

\2((2n + j + \,N)+J2C(k,N) = 

j -=0 Jt = 2 

This follows from Lemma 2(a). (U) 

3 Polylogarithms and related Zeta function 

IDENTITIES 

There exist some interesting identities for the Khinchin constant in terms of 
polylogarithm evaluations. One particularly interesting polylogarithm identity 
is obtained by resolving the integral representation of Corollary 4 in polyloga- 
rithm terms [23]. One may employ the Euler product for si nz/z to write the 
integral as a sum of logarithmic integrals, each in turn expressible in terms of 
polylogarithms. This procedure leads to the series: 


11 A 

log(A'o) log(2) = log 2 (2) + Li 2 (--) + - y ,(-l) n 

L n = 2 

where, the polyiogartthm Li m (z) YllcL =i z k k~ m ■ 

A more direct application of poly logarithms is to invoke the classic Abel identity: 

log(l-a?)log(l— y) - Li 2 (-^— )+Li 2 (— )— Li 2 (a:) — Li 2 (j/)— Li 2 (-r^ r 

1 - y 1 — x (1 - x)(l - y) 

to Lemma 1(c), with x := l/n,y — 1/n to obtain, from a telescoping sum: 

7r 2 1 00 — 1 

log(A' 0 )log(2) = — - - log 2 (2) + J2 Li 2(^T } ' 

n= 2 

An interesting line of analysis starting from this last polylogarithm series is to 
“peel off” parts of the Li 2 function, casting the corrections in closed form. Such 
a procedure gives polylogarithm-based analogues of Theorem 3. For example, 
one can replace the last Li 2 summand above with a more rapidly decaying term: 

r r -i v -1 I L_ 

12 n 2 — 1 n 2 — 1 4 (n 2 — l) 2 

and add back a correction: 


-«(!) + jt!(2)=^ 


59 

64 


where Q is the zeta-like function 
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°o 

n(m) := V ■ ; , V . 

A careful Eulerian partial fraction decomposition (as detailed in [5]) produces 
a closed form for all integral m > 0 


Q(m) 


Lm/2J /o L 


C(2*) 


£) i 

22m-fl 2 


Moreover, 

oo 


n = 1 

and for m = 2,3, • • •, one may show inductively using another partial fraction 
argument 



{C(2n)-i} = n(i) + n(2), 


E ( 2 ) K( 2n ) - U = Q (™) + fi ( m + 1) - (C(2m) - 1}, 

n=l ' ‘ 

from which we may easily obtain a closed form for 

OO OO 

E n*){C(2t) ~ 1} = E^K(2*. !) 

^=1 Jt=l 

for any polynomial V . Thus, 


£(„*)«“■■) = J, 

k = 1 

£(0« 2M) = i? + 5« 2 t 

E( 2 )c( 2 *, 1 ) = || - C(4) + |C(2), 

£(> k ' i) = s~ c(6) n c(4) -^ c(2) - 

E( 4 )c( 2^1) = ^ - C(8) - ^C(4) + j|gC(2), 

and 


*=i 


e :«*.» = 


2069 , . 1 1 7 

2048 - < (l°) + 32 <(6 ) + I 28 C(4 )- 256« 2 ) 


2048 
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with direct analogues for C(2fc, m) for m > 1. 

To complete this section we observe that for m = 1, 2, • ■ ■ 


E 


C(2n) - 1 
n m 


= ^Li m (l/n 2 ). 

n — 2 


Here, unlike Lemma 2(b) in which m = 1, the right-hand sum does not appear 
to telescope. This identity is easy to verify. 


4 Explicit continued fractions 

It is remarkable that, even though a random fraction’s limiting geometric mean 
exists and furthermore equals the Khinchin constant with probability one, not 
a single explicit continued fraction has been demonstrated to have geometric 
mean Aq. Likewise for any negative integer p, not a single explicit fraction has 
been shown to have Holder mean equal to K p . In any event, it is worthwhile to 
mention some explicit fractions with respect to this theoretical impasse. 

The continued fraction for e is 


e=[2, 1,2, 1,1,4, 1,1, 6,1, 1,8,...]. 

The elements are eventually comprised of a meshing of two arithmetic progres- 
sions, one of which has zero common difference while the other has difference 
two and diverges. Thus the meshing has diverging geometric mean. Thus e 
does not possess geometric mean A' 0 . The harmonic mean for e does exist, but 
equals 3/2 which is not A _ i . It turns out that any fraction with elements lying 
in a single arithmetic progression can be evaluated in terms of special functions. 
Explicitly, for any positive integers a,d 


, a T cl y a T 2 d } <2 T 3d, 



where I u is the modified Bessel function of order v. These arithmetic progression 
fractions are certainly interesting, and not beyond deep analysis. It was known, 
for example, to C. L Siegel that these fractions are transcendental [17]. But 
each such fraction has diverging geometric mean and indeed diverging Holder 
means. Note that the means are monotone non- decreasing in p and so a fraction 
with liminf of its elements infinite has infinite means. 


Another example of interest is i r, whose continued fraction expansion is 


* = [3, 7, 15, 1,292, 1, 1,1, 2, 1,3,...] 

The continued fraction elements do not appear to follow any pattern and are 
widely suspected to be in some sense random. Based on the first 17 million 



10 


Bailey, Borwein and Crandall 


digits, the geometric mean is 2.686393 and the harmonic mean is 1.745882 [10]. 
These values are reasonably close to K 0 and AT i , but of course no conclusion 
can be drawn beyond this. 

It is a well known theorem of Lagrange that the elements of a simple continued 
fraction form an eventually periodic sequence if and only if the fraction is an 
irrational quadratic surd. All Khinchin means K p for p - 0,-1, -2,... then 
exist, and are completely determined by the mean of one period of elements. 
Hence each Khinchin mean of a quadratic surd is an algebraic number. Clearly, 
for any algebraic number c = a l ! b formed from integers a, b , one can write down 
a quadratic surd having geometric mean c. Along these lines, it is not hard to 
show that if there exists an integer m > 2 such that 

I°g(A -o/m) 
log(2/m) 

is rational, then there exists a quadratic surd with geometric mean A' 0 . Thus 
the issue of transcendence for Ao and related numbers is an interesting one, 
and one we return to in the next section. It is also of interest that Liouville 
transcendentals form a residual null set, and thus comprise a ’’numerous” set 
whose members typically do not respect the A'o limit. 

If one were in possession of K 0 to arbitrary accuracy, one could of course con- 
struct a fraction having geometric mean I< 0 by appending a “2” (respectively, 
“3”) to the element list whenever the current geometric mean were above (be- 
low) A 0 . There seems to be no way to determine a priori the value of, say, the 
rc-th element. Thus such a constructed fraction is not explicit. 

However, we propose here, without recourse to the value of AT 0 , a continued 
fraction, possessed of explicitly defined elements a n} that should have geometric 
mean Ao- Our construction is based on a deterministic stochastic sampling of 
the Gauss-Kuz’min density, and proceeds as follows. First, for non-negative 
integer n define the van der Corput discrepancy sequence associated with n to 
be the base-2 number 


d(n) — O.60&1&2... 

where the are the binary bits of n> with 60 being least significant. As n 
runs through positive integers, the sequence of d{n) is confined to (0,1) and has 
appealing pseudorandom properties. The construction of the number we shall 
call Z 2 then starts with a 0 := 0, and loops as follows: 

For n = 1 to 00, set a n := [l/(2^ n l — I)J 

The continued fraction elements a n thus determined start out: 
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while the numerical value of Z 2 is approximately 

Z 2 = 0.46107049595671951935414986933669968767808281325747 .... 

On the basis of antic ipated statistical features of our construction we are moved 
to posit: 

Conjecture. The geometric mean of the number Z 2 is in fact the Khinchin 
constant Kq. Furthermore, every p-th Holder mean of Z 2 for p = — 1, -2, ... is 
the respective Khinchin mean K p . 

There are some features of the Z 2 construction that lend credence to our con- 
jecture. On the basis of the known density properties of the discrepancy set, it 
can be shown that the elements of Z 2 are unbounded, and that every possible 
integer element value a n = k > 0 is attained infinitely often. 

With regards to the above conjecture, S. Ploufie [14] has reported a computation 
of the geometric and harmonic means through 5206016 continued fraction ele- 
ments of Z 2 . His results are 2.6854823207 and 1.7454074435, respectively, which 
are remarkably close* to the expected theoretical values. These results raise a 
deeper question: what is the rate of convergence of these empirical means to 
their limiting values' 7 The authors have also been informed by T. Wieting that 
he has an unpublished proof of the basic conjecture, i.e. that the limiting Holder 
means of Z 2 exist for p — 0, —1, —2,... and furthermore equal the corresponding 
Khinchin means [20]. 

In the numerical evaluations reported above and in the next section, the follow- 
ing definitions of empirical Holder means were used: 

/ N \ l f N 

H ° = (n-J 

/ n ytp 

H > ■■= p<° 

Here H p is the quantity that should, as N — + 00 , converge to K p with probability 
one. 

5 Computation of Khinchin Means 

The authors have explicitly computed K 0 and A'_i to more than 7350 decimal 
digit accuracy. These computations were performed with the aid of the MPFUN 
multiprecision softwa-e [2, 3], which was found to be significantly faster for our 
purposes than other available multiprecision facilities. One utilizes this software 
by writing ordinary Fortran-90 code, with multiprecision variables declared to 
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be of type mp.integer, mp_real or mp.complex. In the computations described 
below, the level of precision was sufficiently high that the “advanced” routines of 
the Fortran-90 MPFUN library were employed. These routines employ special 
algorithms, including fast Fourier transform (FFT) multiplication, which are 
efficient for extra-high levels of precision. 

Aq was computed using the formula given above in Theorem 3, with the free 
integer parameter N = 100, and with N = 120 as a check. The implementation 
of this formula was straightforward except for the computation of the Riemann 
zeta function. To obtain 7350 digit accuracy in the final result, 2048 terms of 
the indicated series were evaluated, which requires {C(2fc), 0 < k < 2048} to 
be computed. One approach to compute these zeta function values is to apply 
formulas due to P. Borwein [6]. These formulas are very efficient for computing 
one or a few zeta function values, but when many values are required as in this 
case, another approach was found to be more efficient. This method is based 
on an observation that has previously been used in numerical approaches to 
Fermat’s “Last Theorem” [7, 8]; namely, 

-2 00 

coth(ra) = — y^((2fc)(-l)*:K 2fc 

XX k= 0 

= cosh(7r.r)/ sinh(^x) 

_ 1 1 - f - (7tx) 7 /2\ + (irx) 4 /4\ + (7tx) 6 /6[ H 

7 tx 1 + (7tx) 2 /3! + (7ra:) 4 /5! + (7rx) 6 /7! + ■ • • ‘ 

Let N(x) and D(x) be the numerator and denominator polynomials obtained 
by truncating these two series to n terms. Then the approximate reciprocal 
Q(x) of D(x) can be obtained by applying the Newton iteration 

Qk+i(x) := Qk(x) + [1 - D(x)Q k (x)]Q h (x). 

Once Q(x) has been computed to sufficient accuracy, the quotient polynomial is 
simply the product N(x)Q(x). The required values C(2k) can then be obtained 
from the coefficients of this polynomial. 

Computation time for the Newton iteration procedure can be reduced by start- 
ing with a modest polynomial length and precision level, iterating to con- 
vergence, doubling each, etc., until the final length and precision targets are 
achieved. Computation time can be further economized by performing the two 
polynomial multiplications indicated in the above formula using a FFT-based 
convolution scheme. In our implementation, FFTs were actually performed at 
two levels of this computation: (i) to multiply pairs of polynomials, where the 
data elements to be transformed are the multiprecision polynomial coefficients, 
and (ii) to multiply pairs of multiprecision numbers, where the data elements 
to be transformed are integers representing successive sections of the binary 
representations of the two multiprecision numbers. 
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K- 1 was computed by applying the formula in Corollary 7. Again, the challenge 
here is to pre-compute values of the Riemann zeta function for integer values. 
But in this case both odd and even values are required. The odd values can be 
economically computed by applying the following two formulas, the first given 
by Ramanujan, but simplified slightly, the second derived by differentiating a 
companion identity of Ramanujan ([4]): 


C( 4 yv + 3) =■ -2^ 


1 


£4N+3£ eX p(2fc7|-) — 1) 


k - 1 

2/V+2 

— x(2x) 4Ar+2 £ (-1) 

fc= 0 


k &2k- B^-\-A-2k 
(2A:)!(4A r -{- 4 — 2k)\ 




(27 rk + 2A r ) exp(27r&) — 2N 


k = l 


k* N +'(exp(2kir)~ l) 2 


_ ST'S nfc B 7k BiN+2-2k 

<>N ^ } ^ (2fc- 1)!(4JV+ 2- 2*)!' 


k = 1 


Here B 2 k is as always the 2fc-th Bernoulli number. 

Alternatively the formulas can be written in terms of the even zetas as 


C(4iV + 3) = 


_2 V - 

k 4N+3 (exp(2kr) - I) 

+ 1< (4 ~ 2 +7) C(4^ + 4) - ]T 2C(4*K(4AT + 4 - 4*)} 

k = 1 


1 C ° 

C(4« + 1 ) = --•£ 


(27rk + 2N)exp(27rk) — 2 N 


N ^ k 4N + i (exp(2k 7 r)- l) 2 


k~l 


2 N 


+ 2 - 2*) + (2 N + 1 )C(4N + 2)}. 


k = 1 


These two formulas are not very economical for computing a single odd value 
or just a few odd values of £(k) — again, the formulas in [6] are more efficient 
for such purposes. But these Ramanujan formulas are quite efficient when a 
large number of odd zetas are required. Note that the infinite series in the two 
formulas can be inexpensively evaluated for many N simultaneously, since the 
expensive parts of these expressions do not involve N . Further, the evaluation 
of the infinite series can be cut off once terms for a given N are smaller than 
the “epsilon” of the numeric precision level being used. Happily, convergence 
here is fairly rapid for large N. 

At first glance, the latter summations in these two formulas may appear quite 
expensive to evaluate. But note that each is merely the polynomial product of 
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two vectors consisting principally of even zeta values. Thus both sets of sum- 
mation results can be computed using multiprecision FFT-based convolutions. 

Computation of K 0 to 7350 digit precision required 2.5 hours on an IBM RS6000 / 
590 workstation, and computation of K_i required some 12 hours. Excerpts of 
the resulting decimal expansions for each are included in the appendix. The 
complete expansions are available from the authors. 

One intriguing question is whether the continued fraction expansions of the 
Khinchin constants themselves satisfy the geometric and harmonic mean condi- 
tions of their definitions. We found that the geometric and harmonic means for 
our value of Ao were 2.663660 and 1.746398, respectively. The geometric and 
harmonic means for our value of K_ l were 2.723115 and 1.746965, respectively. 
The issue of where to terminate the list of continued fraction elements is an in- 
teresting one. We employed a simple criterion: if x is known numerically, to D 
decimals to the right of the decimal point, generate continued fraction elements 
for x until a convergent p/q has 2 q 2 >10^. The reason for choosing this simple 
criterion is the theorem that at least one of any two successive convergents must 
satisfy 


and conversely, any reduced ratio p/q satisfying this inequality must be a con- 
vergent of x ([11]). 

To give statistical perspective to these results, we computed these same empiri- 
cal means for 100 pseudorandom multiprecision numbers of the same precision. 
The average and standard deviation of their geometric means were 2.683740 
and 0.030124, respectively. The same statistics for their harmonic means were 
1.745309 and 0.011148, respectively. Note that these two averages are in good 
agreement with the theoretical values K 0 and A'_i. In any event, it appears 
that the empirical geometric and harmonic means for K 0 and A'_i are within 
reasonable statistical limits of the expected theoretical values. 

A question implicitly asked in the previous section is whether A'o or A'_i is al- 
gebraic. This question can be numerically explored by means of integer relation 
algorithms. A vector of real numbers (xi, £ 2 , x n ) is said to possess an 

integer relation if there exist integers a k such that aiXi + a 2 £ 2 H ha n x n = 0. 

It can easily be seen that a real number a is algebraic of degree n — 1 if and 
only if the vector (1, a, a 2 , * * • , aT -1 ) possesses an integer relation. Even if a 
is not algebraic, integer relation algorithms produce bounds that allow one to 
exclude relations within a region. 

We employed the “PSLQ” algorithm developed by Ferguson and one of the 
authors, a simplified version of which is given in [1]. This algorithm, when 
applied to power vectors generated from our computed values of A r 0 and A\_i, 
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found no relations for either. On the contrary, we obtained the following result: 
neither K$ or K_i satisfies a polynomial of the form 

0 — ao 4- ai<* 4 a 20r 2 4 o,$o? 4 ■ • • 4 as oo 50 

with integer coefficients a * of absolute value 10 70 or less. 

In a second experiment, we explored the possibility that K 0 or K_ v is given by 
a multiplicative formula involving powers of primes and some well-known math- 
ematical constants. To that end, let denote the &-th prime. We established, 
using PSLQ, that neither Ao or A'_i satisfies a relation of the form 

15 

0 = q 0 log or 4 g k log p k 

k= 1 

4 a i6 log 7 r 4 a 17 loge 4 a. 18 log 7 4 a 19 log ((3) 4 a 2 o log log 2 

with integer coefficients a& of absolute value 10“^ or less. By exponentiating 
this expression, it follows that neither Ao or A_i satisfies a corresponding mul- 
tiplicative formula with exponents of absolute value 10 20 or less. 

Acknowledgment. Thanks are due to Robert Corless, Greg Fee, Thomas 
Wieting, Simon PloufFe, and Joe Buhler for many helpful discussions. 
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Appendix: The Khinchin Constant A' 0 to 7,350 Digits 

2 . 

68545200106530644530971483548179569382038229399446 

29530511523455572188595371520028011411749318476979 

95153465905288090082897677716410963051792633483259 

66838185231542133211949962603932852204481940961806 

86641664289308477880620360737053501033672633577289 

04990427070272345170262523702354581068631850103237 

46558037750264425248528694682341899491573066189872 

07994137235500057935736698933950879021244642075289 

74145914769301844905060179349938522547040420337798 

56398310157090222339100002207725096513324604444391 


36909874406573435125594396103980583983765664559601 

The Khinchin Harmonic Mean A'_i to 7,350 Digits 
1 . 

74540566240734686349459630968366106729493661877798 

42565950137735160785752208734256520578864567832424 

20977343982577985596531102601834294460206578713176 

15026238960612981165718728271638949622593992929776 

06160830078357479801549029312671643067241248453710 

96077711207484391474195803753220015690822609477078 

44894635568203493582068440202422591615018316479048 

29229656977733143662210991806388842581650599997697 

61391683577259217628635718712601565066754443340174 

00283376465305136584406098398017126202832041200630 


78553128249666473680304034761497467330708479436280 


Khinchin Means A p for Various Negative p to 50 Digits 
P K P 

-2 1 . 450340328495630406052983076680697881408299979606904 . 
-3 1 -313507078687985766717339447072786828158129861484792. 
-4 1-236961809423730052626227244453422567420241131648937. 
-5 1 . 189003926465513164062363732771403397386092512639671 
-6 1 . 156552374421514423152605998743410046840213070718761 
-7 1 . 133323363950865794910289694908868363599098282411797 
-8 1 . 115964408978716690619156419345349696769491182230400 
-9 1-102543136670728013836093402522568351022221284149318 
-10 1 . 091877041209612678276110979477638256493272651429656 . 
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